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The deconfinement transition in 3+1 dimensional gluodynamics is studied using the gauge invariant variational method 
introduced by Kogan and Kovner a few years ago. We identify a first order phase transition, characterized by a dis- 
continuous jump in the entropy of the system, resulting in a transparent picture of the mechanism of deconfinement. 
The calculation of the ratio of the transition temperature to the mass of the lightest glueball in the model yields 0.18 in 
complete agreement with the lattice estimate. 



1 Introduction 

In the almost 25 years since the pioneering work of 
Polyakov [Q] and Susskind much effort has been 
devoted to attempts to understand both the basic 
physics and quantitave features of the deconfining 
phase transition of QCD. 

The high temperature phase, i.e. for temperatures 
well above the critical temperature T c , is becoming 
well understood, and is widely believed to resemble a 
plasma of almost free quarks and gluons. However, 
the transition region, T c < T < 2T Cl is very poorly 
understood. This region is the most interesting one 
since it is there that the transition between 'hadronic' 
and 'partonic' degrees of freedom occurs. 

The study of the transition region is a complicated 
and inherently non-perturbative problem which has 
mostly evaded treatment by analytical methods. The 
variational method introduced several years by Kogan 
and Kovner appears well suited to a discussion of 
such a problem. 

Recently, the method was applied to the study of the 
phase transition in SU (N) gluodynamics. As a result 
of these studies [IHE], we identified a first order phase 
transition at a transition temperature, in units of the 
lighest state in the theory, of T c = 0.18, in complete 
agreement with the lattice estimate for SU(3) [Ej. 

We minimize the relevant thermodynamic potential at 
finite temperature, i.e. the Hclmholtz free energy, on 
a set of suitably chosen trial density matrix function- 
als. As with any variational calculation there is an 
element of guesswork involved. However, the choice 
of the trial states is dictated by rather general prin- 
ciples. As in the quantum mechanical case, the trial 
state ought to be simultaneously rich enough to span 



the interesting physics and to lead to calculable re- 
sults. This last requirement is particularly stringent, 
since only functional integrals of the gaussian form can 
be explicitly evaluated. Specific to the quantum field 
case are the requirements of gauge invariance and of 
correct UV behaviour. To understand this last state- 
ment, it suffices to mention first, that the free energy 
is dominated by UV modes while we are interested in 
the IR physics and, second, that in an asymptotically 
free theory such as gluodynamics, low and high mo- 
mentum modes are, in general, non-separable. The 
method of [ |S] to be discussed below accommo- 
dates all these requirements, leading to what is, to the 
author's best knowledge, the only first principle ana- 
lytical analysis of the deconfinement phase transition 
in gluodynamics. 



2 The variational Ansatz 

We consider density matrices which, in the field basis, 
have Gaussian matrix elements and where gauge in- 
variance is explicitly imposed by projection onto the 
gauge-invariant sector of the Hilbert space 



p[A,A 







j DU exp | 









A1(x)G- la \x,y)A b Ay) 



+ A i Ua (x)Gr j lab (x,y)A J ub (y) 
-2A?(x)H$(x,y)A' ub (y)]\, (1) 



*In memory of Ian Kogan: friend, mentor and colleague. 



where f xy = f d 3 xd 3 y, DU is the SU(N) group- 
invariant measure, and under an SU (N) gauge trans- 
formation U 

Af(x) - A?*(x) = S ab (x)A b (x) + \?(x) , (2) 
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with 



S ab = -tr(T°UV[/) , 

A? = -tx(r a U^diU) , 
9 



(3) 
(4) 



and \- form an N x N Hermitian representation 

of SU(N): [y, 2 ^] = if abc ^ with normalization 
tr (r Q r b ) = 2<5 a, > 

We take the variational functions diagonal in both 
colour and Lorentz indices, and translationally invari- 
ant 



G 



-lab 



(x,y) = S ab 6 lJ G- 1 (x-y) 
H? b {x lV ) = 5 ab 5^H(x-y). 



(5) 



Further, we split the momenta into high and low 
modes with k M and restrict the kernels G _1 and 
H to the one parameter momentum space forms 



G- 1 (k) = 



M, k < M 
k, k > M 



H(k) i H > k < M 

H[k >- < \ 0, k>M 



(0) 



The logic behind this choice of ansatz is the following. 
At finite temperature we expect H(k) to be roughly 
proportional to the Bolzmann factor exp{ — E(k)(3}. In 
our ansatz, the role of one particle energy is played by 
the variational function G _1 (fc) and its form is mo- 
tivated by the propagator of a massive scalar field, 
i.e. (k 2 + M 2 ) 1 / 2 . We will be interested only in tem- 
peratures close to the phase transition, and those we 
anticipate to be small, T c < M . For those tempera- 
tures one particle modes with momenta k > M are not 
populated, and we thus put H(k) = 0. For k < M the 
Bolzmann factor is non- vanishing, but small. Further, 
it depends only very weakly on the value of the mo- 
mentum. With the above restrictions on the kernels, 
only two variational parameters, M and H, remain. 

Importantly, the density matrix functional in eq. (fll) 
describes, for H = 0, a pure state p = >< 
where ^[A] are Gaussian wave functionals. For H ^ 0, 
eq. Q describes a mixed state with \H\ proportional 
to the entropy of the trial density matrix. 

The expectation value of a gauge invariant operator in 
the variational state eq. is then given by 

(0) A ,u = Z- l Tr{pO) 



VUVAO(A,A')- 



exp 



AG A 



A^G-'A' 1 



2AHA' L 



(7) 



where Z is the normalization of the trial density ma- 
trix p, i.e. 

Z = Tip 

VUVA- 



exp 



AG^A + A U G- 1 A U 



2AHA 



u 



(8) 



To evaluate the above expressions we first perform, for 
fixed U(x), the gaussian integration over the vector 
potential A. For Z we get, in leading order in H, 

z= fvu c ^l- l -x{^ + ^(s + s T ))x 



-HGtr(S + S')}. (9) 



We now integrate out the high momentum, k 2 > M 2 , 
modes of U perturbatively to one-loop order. This ef- 
fects a renormalization group transformation on the 
low modes, replacing the bare coupling g 2 by the run- 
ning coupling g 2 (M) [EJ[7]. To one- loop accuracy, 
the coupling g 2 (M) runs identically to the Yang-Mills 
coupling 

3 The effective a-model 

The normalization Z can be then interpreted as the 
generating functional 



Z = Tip = J VUe- s{u) 



(10) 



for an effective non-linear tr-model for the low momen- 
tum modes (k 2 < M 2 ) defined by the action 

^tr(dUdU^) 

Uu^dU - dU j U)(dUU^ - UdU r ) 



S(U) 



H 



8g 2 
1 

i^ 2 



tr 



HMhrUhrU . 



(11) 



where U independent pieces have been dropped. 



The matrix U plays the same role as Polyakov's loop 
P at finite temperature — the functional integration 
over U projects out the physical subspace of the large 
Hilbert space on which the Hamiltonian of gluodynam- 
ics is defined. 

This (T-model has a phase transition [j5| at the critical 
point (for Kqcd = 150Mev, N = 3 and with the one- 
loop Yang-Mills f3 function) 



Mr = A 



QCD 



e 11 



.86A 



QCD 



1.33Gev. 



(12) 
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For M < M c , the a- model is in a disordered, 
SU(N)l®SU (N)r symmetric, phase with massive ex- 
citations and where (U) — 0. Since U is the Polyakov 
loop, this corresponds to a confined state. When 
M > M c , the cr-model is in an ordered, SU(N)y 
symmetric, phase with massless Goldstone bosons for 
which (U) ^ 0, corresponding to a deconfined state. 
With this analysis we have established a correspon- 
dence between the cr-model phase transition and the 
deconfinement transition in SU(N) gluodynamics. 

In fact, this correspondence can be argued in rather 
general terms. Instead of restricting ab initio the den- 
sity matrix to the form eq. imagine that we take 
some arbitrary gauge-invariant density matrix ansatz 
that depends on the A fields and is integrated over the 
U fields. We allow this new ansatz (and whatever ker- 
nels it may contain) to remain arbitrary until we have 
no choice but to restrict it. Then we integrate out the 
A fields to obtain a partition function of U fields with 
respect to some action. 

Next we introduce a separation of momenta into high 
and low modes with k ^ M and integrate out the high 
mode U fields as before. This effects a renormalisation 
group transformation on the low modes, replacing the 
bare coupling g 2 — which is not arbitrary, since it is 
defined by the gauge transformations eq. J2J — by the 
running coupling g 2 {M). Now provided our ansatz 
is sufficiently close to the correct density matrix for 
SU(N), the theory will be asymptotically free. We 
are thus left with an action for the low modes which 
is, once again, a complicated cr-model, with a renor- 
malised coupling g 2 (M) which we expect to be small 
provided M is large, and vice versa. 

Now consider this model as a statistical mechanical 
model at 'temperature' g 2 (M). We make the plausi- 
ble assumption that this cr-model will, as M is varied, 
undergo a symmetry-breaking transition at 'tempera- 
ture' g 2 (M c ) fro m a 'thermally disordered' (symmet- 
ric) phase at large g 2 (M c ) to an ordered phase at small 
g 2 {M c ). Further, it is clear — since the Polyakov loop 
(U) is zero in the former phase and non-zero in the lat- 
ter — that this CT-model phase transition corresponds 
directly to the deconfinement transition in the SU (N) 
theory. 

On review, it is clear that our only assumptions are 
that the ansatz is sufficiently close to SU(N) and 
that the low mode cr-model undergoes a symmetry- 
breaking phase transition. In particular, let the 
ansatz, which is arbitrary and need not be Gaussian, 
be the correct density matrix for SU(N). The first as- 
sumption is certainly true. If the second assumption 
is also true, then we have constructed an exact ar- 
gument that the deconfinement transition in SU (N) 
corresponds to the phase transition in the low mode 



cr-model. 

Thus, in order to study deconfinement in SU(N), our 
aim should be to model the physics of each cr-model 
phase as accurately as possible and calculate the tran- 
sition scale M c . We then calculate the free energy of 
SU(N) in each phase, including any possible contri- 
bution from the high modes, at temperature T and 
extract the minimal free energy. The deconfinement 
transition occurs at the temperature for which the 
free energies calculated in the ordered and disordered 
phases of the low mode cr-model coincide. 

Although we will take eq. Q as the ansatz for the 
density matrix, we shall keep the kernels G _1 and H 
arbitrary until we have no choice but to restrict them. 

4 Calculation of the free energy 

The Helmholtz free energy F of the density matrix p 
is given by 

F = (H)-T(S), (13) 
where H is the standard Yang-Mills Hamiltonian 

1 



H 



ii.2 



(14) 



with 



'8A?(x) ' 
1 



B?{x) = -e ijk {d j Al{x)-d k A a j {x) 



(15) 



+ gf ab °A b j (x)A%(x)} l 

S is the entropy, and T is the temperature. 
Thus 

F = i(Tr(£ 2 e) + Tr(B 2 c,)) + T- Trying). (16) 

In the disordered phase, no progress seems possible 
without restricting the arbitrary kernels. Following [ 
0|, we adopt the forms eq. iJBJ. For small H, we con- 
sider only the first non-trivial order in H, that is a 
term of o(H In H ) in the entropy. This term can be 
written as a product of left SU{N) and right SU(N) 
currents and does, therefore, vanish in the disordered, 
SU(N) L ® SU(N) R symmetric, phase The re- 

maining contribution to the free energy, the average 
of the Hamiltonian, is evaluated in the mean field ap- 
proximation [ Ej . The free energy is minimized for 
M = M c ~ 1.33GeV 



Fdi 



N 2 M A C 

307T 2 



(17) 
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The simplest option to evaluate the free energy in the 
disordered phase is to use perturbation theory. Per- 
turbation is certainly appropriate for large enough val- 
ues of M, where the expectation value of the U field 
is close to unity. From numerical studies [ it is 
known that the transition occurs when the expecta- 
tion value of U is greater than .5. We can thus expect 
perturbation theory to be qualitatively reliable down 
to the transition point. In the leading order perturba- 
tion theory approximation to the ordered phase of the 
cr-model, however, minimisation with respect to arbi- 
trary kernels G~ l and H for both high and low modes 
is possible. Further, the analysis can be carried out to 
all orders in the thermal disorder kernel H . 

In this approximation, the U matrices can be parame- 
terised in the standard exponential form and expanded 
in the coupling g 



= 1 



■HHP y 



(18) 



Hence at leading order one can take 
(7-1, 

diU ~igdi<p a ^-. (19) 

Thus, the gauge transformations @ reduce to 

Af -» Af - dnp a (20) 
and the Hamiltonian (|14fl reduces to 



H=\[Ef + (e ijk d j A%) 2 ] 



(21) 



These last two equations describe the theory 
U(1) N _1 : in the leading order of a- model perturba- 
tion theory, the SU(N) Yang-Mills theory reduces to 
the U(l) _1 free theory. Moreover, the density ma- 
trix eq. becomes Gaussian again, because the gauge 
transformations are linear. One has 



p[A,A] = 



AG' 1 A 



+ (A' -d<p)G~\A! -dip) 

- 2AH{A' - dtp) 



(22) 



Now the theory of N 2 — l U{1) free fields in 3+1 dimen- 
sions is completely tractable; the variational analysis 
for the U(l) theory (with Gaussian ansatz H220 ) was 
discussed in [ lll)j . The free energy in momentum space 



in terms of the arbitrary kernels G 1 and H is 

n 2 - 1 r d 3 P 



F =- 



2 J (2tt) 3 
4T( log 



G- X {1 + GH) + P 2 G{1 - GH)- 1 
GH 



- log 



(1 - {GH) 2 ) 1 / 2 - (1 - GH)\ 
1 - (1 - {GH) 2 ) 1 / 2 



GH 

1 - {1 - {GH) 2 ) 1 / 2 
~{l-{GH) 2 y/ 2 -(1-GH). 

The kernels which minimise the free energy are 
H = 2p 



(23) 




(24) 



and the minimal value of the free energy at tempera- 
ture T is 

n 2 N 2 T 4 



F = — 



45 



(25) 



So we see that the free energy of SU{N) is minimised 
with M — M c in the disordered phase of the cr-model 
for temperatures from zero up to a temperature T c 
where 



F = 



N 2 M* 



30tt 2 

which in turn implies 



ir 2 N 2 T* 
45 



7; = [^V /4 ^£~ 470MeV. 



(26) 



(27) 



5 Conclusions 

We find that the deconfinement phase transition is 
strongly first order with a transition temperature of 
470MeV. The transition is due to an abrupt jump in 
the entropy. In the disordered phase it is zero since 
glueballs are heavy and, therefore, their contribution 
is suppressed by the Boltzmann factor exp(— M g /T). 
In the ordered phase, the entropy is non- vanishing and 
proportional to the number of coloured states. 

The major uncertainty in the calculation is due to the 
fact that the scale M c is very sensitive to the mean 
field approximation. It is, therefore, only meaningful 
to compare the results of our approach with lattice 
estimates for M c independent quantities. For one such 
quantity, the transition temperature in units of the 
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mass of the lightest glueball in the model, i.e. 2M C [ 
EJ, we find 



i-(|)V*~ .18. (28) 



2M C 2?r v 2 ' 
which is in agreement with the lattice result. 
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